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Abstract 

We investigated the double-exchange spin ice model on a kagome lattice by Monte Carlo simulation in order to 
study a mechanism for lifting the degeneracy in frustrated electron systems. We show specific heat and vector spin 
chirality data on a finite lattice. Specific heat has a double-peak structure: A broad peak and a sharp peak are at 
k^T/t ~ 0.15 and 0.015, respectively, where t is the transfer integral of electrons. The broad peak corresponds to 
a crossover to the spin ice-like state, on the other hand, the sharp one a transition to a dodecagonal spin cluster 
(dodecamer) state. We discuss the interplay between the formation of the dodecamer state and the lifting of the 
macroscopic degeneracy. 
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Macroscopic degeneracy is one of the well-known 
properties in frustrated systems. However, such a de- 
generacy should be lifted at T = in realistic systems 
because of the third law of thermodynamics. Consider- 
ing these circumstances, it is significant to investigate 
mechanisms to lift the degeneracy. In particular, a pe- 
culiar mechanism is expected in electron systems with 
the frustration: The motion of electrons plays an im- 
portant role for lifting macroscopic degeneracy, which 
may lead to some kinds of exotic phases. 

In order to see such a mechanism, we investigate a 
double-exchange spin ice (DESI) model on a kagome 
lattice by a Monte-Carlo (MC) simulation [1], where 
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the kagome lattice consists of corner-shared up- and 
down-triangles. The DESI model is the Anderson- 
Hasegawa model [2] with a following condition: Lo- 
calized spins have uniaxial anisotropies that they are 
forced to point either inward (in-spin) or outward 
(out-spin) for an up-triangle. In this case, localized 
spins are coplanar. The uniaxial anisotropies and fer- 
romagnetic interaction between the nearest-neighbor 
(n.n.) spins due to the Anderson-Hasegawa mecha- 
nism cause the frustration in the model for the same 
reason seen in the spin ice systems [3] . 

From our previous MC calculations, following re- 
sults have been obtained [1]: At sufficiently low tem- 
peratures, a dodecagonal localized spin cluster, "dode- 
camer" [see Fig. 1] , is realized in a wide doping region 
n ~ 1/3 ~ 1/2, where n is the number of particles per 
site. The dodecamer order is driven by both the kinetic 
energy gain and the frustration. 

In this paper, specific heat and vector spin chirality 
data of the system are shown, in wide temperature 
range in order to investigate a relation between the 
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Fig. 1. Dodccamcr formation in the kagome lattice: The 
shadowed dodccagonal spin cluster is a dodccamcr, where 
in-(out-)spins for up-triangles are represented by open (closed) 
circles. There are two types of bond in the DESI system. One is 
a bond which connects an in-spin and an out-spin (solid line), 
and the other is that which connects the n.n. pairs of the same 
spins (dashed line). The electron transfer integral has larger 
value on the former bond than the latter. The dodecamer con- 
sists of twelve spins connected by solid bonds along which elec- 
trons can move easily. The dodccamcr order is a bond order 
rather than a spin order. The dodecamer state means transla- 
tional symmetry breaking. 

entropy release and the spin structure. According to 
the results obtained here, it is found that the system 
undergoes spin ice- like states as in the spin ice systems 
[3] at intermediate temperatures and the dodecamer 
state at further low temperatures. 

We have performed MC calculations [4] to study the 
thermodynamics of the system at finite temperatures. 
We typically run 100,000 MC steps for measurement 
after 10,000 thermalization steps. 
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Fig. 2. Temperature dependence of the specific heat and a 
vector spin chirality (v z ) for 6x6 unit cells (— 108 sites) and 
the chemical potential — 0. Error bars are within the sizes of 
the symbols. Specific heat shows a double-peak structure: One 
is at k B T/t ~ 0.015 and the other k B T/t ~ 0.15. (v z ) changes 
rapidly at around the broad peak of the specific heat data. 
Furthermore, the sharp peak exists in the (v z ) ~ 1 region. 

In order to see how the entropy is released in the 
system, specific heat was investigated, and data are 
shown in Fig. 2 for 6 x 6 unit cells (= 108 sites) with 



the chemical potential /i = 0. We can see a double-peak 
structure: A broad peak is seen at ksT/t ~ 0.15 and a 
sharp peak is at kBT/t ~ 0.015, where t is the transfer 
integral of electrons. 

To clarify the origin of the broad peak, we check a 
spin configuration of the system. A vector spin chirality 
is defined for triangles in the kagome lattice which has 
the form, v = --^g (s 4 x Sj + Sj x&l&x Si) , 
where i, j and k represents three sites of a trian- 
gle. Here, we set \Si\ = 1. We only consider the 
^-component of v, because localized spins are copla- 
nar in the system. Each triangle is allowed to take 
"three-out" (v z = —3), "three-in" (v z — —3), "two- 
in one-out" (v z — 1) or "one-in two-out" (v z = 1) 
spin structures. At sufficiently high temperatures, all 
of these patterns are realized equally, i.e., (v z ) = 0, 
since both one-in two-out and two-in one-out have a 
three-fold degeneracy. On the other hand, as the tem- 
perature is lowered, three-out and three-in states are 
excluded, because they have higher energy than the 
others due to the n.n. ferromagnetic interaction. Then, 
the system has (v z ) — > 1. This situation is similar to 
the spin ice systems [3]. Thus, we call such states an 
"ice state", hereafter. 

Temperature dependence of (v z ) are shown in Fig. 2. 
(v z ) grows as the temperature is lowered, and almost 
saturates below the broad peak. This means that the 
system becomes the ice state at the broad peak tem- 
perature. Thus, we conclude that the broad peak of the 
specific heat data corresponds to a crossover to the ice 
state. Note that the ice state still has the macroscopic 
degeneracy due to the frustration. 

Next, let us consider the sharp peak in the specific 
heat data. We investigated a dodecamer structure fac- 
tor D q in Ref [1]. The largest peak of D q is defined 
as D q , which shows a rapid change at ksT/t ~ 0.015. 
Temperature dependence of D q indicates that there 
is an inflection point in this temperature range. The 
sharp peak temperature is almost consistent with that 
of the inflection point. Thus, it is considered that the 
sharp peak corresponds to a phase transition from the 
ice state to the dodecamer state, which is smeared out 
by the finite-size effect. 

In conclusion, the process for lifting the degeneracy 
in the system is as follows: At intermediate tempera- 
tures, the n.n. ferromagnetic interaction leads the sys- 
tem to the ice state, which still has the macroscopic 
degeneracy. According to this result, the short-range 
correlation is considered to be insufficient to lift the de- 
generacy, often seen in the classical spin systems. Fur- 
thermore, at sufficiently low temperature region, the 
effective long-range interaction is relevant in order to 
gain the kinetic energy, and the degeneracy is lifted by 
the dodecamer formation accompanied with the trans- 
lational symmetry breaking. 
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